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Abstract: Currently, the long-standing issue concerning the size of the gluonic content 
of the r]' and r] mesons remains unsettled. With this in mind we consider the central 
exclusive production (CEP) of ij' , rj meson pairs in the perturbative regime, applying the 
Durham pQCD-based model of CEP and the 'hard exclusive' formalism to evaluate the 
meson production subprocess. We calculate for the first time the relevant leading order 
parton-level processes gg — )• qqgg and gg — )• gggg, where the final-state gg and qq pairs 
form a pseudoscalar flavour-singlet state. We observe that these amplitudes display some 
non-trivial and interesting theoretical properties, and we comment on their origin. Finally, 
we present a phenomenological study, and show that the cross sections for the CEP of rj' , 
T] meson pairs are strongly sensitive to the size of the gluon content of these mesons. The 
observation of these processes could therefore provide important and novel insight into this 
problem. 
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1. Introduction 



Central exclusive production (CEP) processes of the type 

pp{p) ^P + X + p{p) , (1.1) 

can significantly extend the physics programme at high energy hadron colliders. Here 
X represents a system of invariant mass Mx, and the '+' signs denote the presence of 
large rapidity gaps. Such reactions provide a very promising way to investigate both QCD 
dynamics and new physics in hadron collisions, and consequently they have been widely 
discussed in the literature, with recently there being a renewal of interest in the CEP 
process, sec for example [1-3] for reviews and further references. 

It is well known that the r] and r]' mesons, the isoscalar members of the nonet of the 
lightest pseudoscalar mesons, play an important role in the understanding of low energy 
QCD. Knowledge of the quark and gluon components of their wave functions would provide 
important information about various aspects of non-perturbative QCD, see for instance [4] 
and references therein. A clear observation of the presence of a purely gluonic, gg, com- 
ponent in the r]' (and r]) mesons would also confirm that the gluons play an independent 
important role in hadronic spectroscopy. The presence of such a component of the t]' me- 
son, due to the so called gluon anomaly [5-8] , is also related to the old question of why the 
rj' mass is much larger than that of the rj (the well-known '^7(1) problem', see [9]). 

Currently, while different determinations of the rj-ri' mixing parameters are generally 
consistent, the long-standing issue concerning the extraction of the gluon content of the 
rj' (and rj) remains uncertain, in particular due to non-trivial theory assumptions and 
approximations that must be made, as well as the current experimental uncertainties and 
limitations, see for example [10] for a discussion of the theoretical uncertainties (in e.g. the 
decay form factors) present in such extractions and [4] for a review of the experimental 
situation. The results of a detailed analysis of various radiative processes and heavy particle 
decays, see for example [4,11,12], therefore do not currently allow a conclusive confirmation 
or otherwise of a non-gg component in the the r]' and r] mesons, within the experimental 
uncertainties. Moreover, as discussed in [4], it is unlikely that lattice simulations of QCD 
will provide a determination of the gluonic contribution of the r]', r] wave functions in the 
near future. On the other hand, in [13] a study of the 777 and rj'^ transition form factors, 
F(^^/)^((5^), appears to indicate that the two-gluon Fock component of the rj' meson may 
be quite large. 

In this paper we will show that the CEP of pseudoscalar meson pairs {ij't]' , rjij' , r]r]) at 
sufficiently high transverse momenta p± can provide a potentially powerful tool to probe 
the structure of the r]', r] mesons, and is especially well suited to addressing the old problem 
of the value of the gluonic flavour-singlet contribution, as well as clarifying the issue of 77-77' 
mixing. We show that any sizeable gg component of the 77' (and 77) can have a strong effect 
on the CEP cross section, and therefore such an exclusive process represents a sensitive 
probe of this. Moreover, we may expect that in the near future, after analysing the existing 
4 photon candidates with Et > 2.5 GeV and forward rapidity gaps, CDF will collect a large 



- 1 - 



number of ry',?] events and also significantly improve the current limits on vr^vr'^ CEP [14]. 
T] and rj' meson CEP may also be studied within the CMS/TOTEM special low-pileup 
runs with sufficient luminosity [14, 15]. 

In [16] the observation of 43 77 events with [77(7)] < 1.0 and Et{'j) > 2.5 GeV, with 
no other particles detected in — 7.4 < < 7.4 was reported, which corresponds to a cross 
section of a^^ = 2A8toil is^^^)-oii (syst) pb. In [17] the 7r°7r° 4j CEP background was 
calculated for the first time and found to be small (with a{7r^7r^)/a{'y^) ~ 1%), a prediction 
that is in agreement with the CDF measurement, which finds that the contamination caused 
by 7r°7r° CEP is very small (< 15 events, corresponding to a ratio N{Tr'^7T^)/N{jj) < 0.35, 
at 95% C.L.). As well as representing a potential observable, we note that ri{')r]{') CEP, 
via the r/C) — >■ 77 decay may also in principle represent a background to 77 production, 
if this cannot be suppressed experimentally. It is therefore important to calculate the 
predicted ??(')'?(') CEP cross sections, in particular in the presence of a potentially large 
gg flavour-singlet component, which may enhance the corresponding rates. 

We will apply the 'hard exclusive' formalism described in [18] (see also [19]) to the 
production of a meson, M anti-meson, M, pair via the gg MM subprocess, which may 
then be used to calculate the corresponding CEP cross section, for X = MM in (1.1). 
In [17] the gg — t- qqqq parton-lcvcl amplitudes were calculated, where each qq pair forms a 
meson state, and from this the flavour singlet and non-singlct meson pair CEP cross section 
was calculated. However, any non-zero gg component of the flavour singlet ry' and, through 
mixing, 77 mesons will also be accessed via the gg — > ggqq and gg — > gggg subprocesses. We 
will show that these amplitudes, which on the face of it are quite unrelated, in fact display 
a striking similarity, being identical with each other and with the 'ladder-type' gg — >■ qqqq 
amplitude calculated in [17], that only contribute for flavour-singlet mesons, up to overall 
normalization factors. We show how this remarkable result may be explained in the MHV 
framework by the fact that the same external parton orderings contribute in all three cases. 

Theoretical studies of meson pair CEP in fact have a long history, which predates the 
perturbative Durham approach depicted in Fig. 1. Exclusive tttt production, for example, 
mediated by Pomeron-Pomeron fusion, has been a subject of theoretical studies within 
a Regge-pole framework since the 1970s (see, for instance [20-22] for early references 
and [23-25] for more recent ones). However, as discussed in [24,25], at comparatively 
large meson transverse momenta, k±, CEP should be dominated by the perturbative 2- 
gluon exchange mechanism discussed above and shown in Fig. 1. At lower k± a study of 
the transition region between these 'non-perturbative' and 'perturbative' regimes may be 
necessary, as was performed in [25] for the case of tttt CEP. For the {E± < 2.5 GeV, \r]\ < 1) 
meson pair event selection we will consider in this paper, the perturbative contribution was 
found to be dominant, and this is certainly expected to be true for the r]{')rj{') perturbative 
cross sections, which we will show to be enhanced relative to 7r°7r° production. We will 
therefore neglect such a Regge-based non-perturbative contribution throughout this paper. 

The outline of this paper proceeds as follows. In Section 2 wc introduce the CEP 
formalism for the process (1.1). In Section 3 we introduce the 'hard exclusive' formalism 
used to model the gg — t- MM subprocess amplitudes, and calculate for the flrst time the 
gg ggqq and gg — )• gggg amplitudes through which the gg component of the r]' and rj is 
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accessed in the CEP process. In Appendix B we show how these ampUtudes may also be 
calculated using the MHV formalism, in a way which sheds some light on the interesting 
theoretical features which these amplitudes display. This aims to provide some theoretical 
insight into these amplitudes, but can be skipped by the reader who is only interested 
in the phenomenological implications of our analysis. In Section 4 we present numerical 
results for the CEP of rfT]' and r/r/ meson pairs, for a range of different sizes of the gg 
component of the flavour-singlet distribution amplitude. Finally, we conclude in Section 5. 

2. Central exclusive production 

The formalism used to calculate the perturbative CEP cross section is explained in detail 
elsewhere [26] and so we will only review the relevant aspects here (for recent reviews and 
references see [1-3]). The amplitude is described by the diagram shown in Fig. 1, where 
the hard subprocess gg ^ X is initiated by gluon-gluon fusion and the second t-channel 
gluon is needed to screen the colour flow across the rapidity gap intervals. We can write 
the 'bare' amplitude in the factorized form [25,27] 

T = Tr'^ [ — g-TT ^ ^\2'^ _L_ ^ fg{xi,x[,Ql,fj.'^;ti)fg{x2,x'2,Ql,fi'^;t2) , (2.1) 

where the /^'s in (2.1) are the skewed unintegrated gluon densities of the proton: in the 
kinematic region relevant to CEP, they are given in terms of the conventional (integrated) 
densities g{x,Qf). is the 4-momentum transfer squared to proton i and n is the hard 
scale of the process, taken typically to be of the order of the mass of the produced state: 
as in [27], we use /x = Mx/2 in what follows. The t-dependence of the /^'s is isolated in a 
proton form factor, which we take to have the phenomenological form F]\f{t) = exp(6t/2), 
with b = 4GeV^^. The A4 is the colour-averaged, normalised sub-amplitude for the 
gg ^ X process 




Figure 1: The perturbative mechanism for the exclusive process pp ^ p + X + with the eikonal 
and enhanced survival factors shown symbolically. 
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Here a and b are colour indices, Mx is the central object mass, V^j^ represents the gg ^ X 
vertex and qi_^ are the transverse momenta of the incoming gluons, given by 

QijL =Q±- Pi± > 92x = -Q± - P2x , (2-3) 

where Q± is the momentum transferred round the gluon loop and pij_ are the transverse 
momenta of the outgoing protons. Only one transverse momentum scale is taken into 
account in (2.1) by the prescription 

Qi = min{Q_L, |(Q± - Pix)l} > 

Q2 = min{Q^,|(Q^ + P2^)|}. (2.4) 
The longitudinal momentum fractions carried by the gluons satisfy 

where x' is the momentum fraction of the second t-channel gluon. The differential cross 
section at X rapidity yx is then given by 

/a2 \ f -,2 ,2 l^(PlxiP2 



{SU) J d^Plxd^P2x^^^^^^5lk(Plx,P2x) , (2.6) 



where T is given by (2.1) and S^^^ is the 'eikonal' survival factor, the probability of pro- 
ducing no additional particles due to soft proton-proton rescattering. This is calculated 
using a generalisation of the 'two-channel eikonal' model for the elastic pp amplitude (see 
[28] and references therein for details). 

Besides the effect of eikonal screening, 5eik, there is an additional suppression caused by 
the rescatterings of the intermediate partons (inside the unintegrated gluon distribution, 
fg). This effect is described by the so-called enhanced Reggeon diagrams and usually 
denoted as S^^, see Fig. 1. The value of S^^^ depends mainly on the transverse momentum 
of the corresponding partons, that is on the argument Q'^ of fg{x,x',Qf,fi'^;t) in (2.1), 
and depends only weakly on the pj_ of the outgoing protons (which formally enters only 
at NLO). While in [26, 29] the S^^-factoT: was calculated using the formalism of [30], 
here, following [17,25], we use a newer version of the multi-Pomeron model [31] which 
incorporates the continuous dependence on and not only three 'Pomeron components' 
with different 'mean' Qi. We therefore include the Senh factor inside the integral (2.1), 
with (-Sgj^jj) being its average value integrated over Q±. 

If we consider the exact limit of forward outgoing protons, pij^ = 0, then we find that 
after the Q± integration (2.2) reduces to 



Ai,A2 



M « g\,giVij ^ -Ql{V++ + V..) ~ V 6^^^'Wx,x, , (2-7) 



where A^^ 2) are the gluon helicities in the gg rest frame. The only contributing helicity 
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amplitudes are therefore those for which the gg system is in a = state, where the 
z-axis is defined by the direction of motion of the gluons in the gg rest frame, which, 
up to corrections of order ~ q]_/M'^, is ahgned with the beam axis. In general, the 
outgoing protons can pick up a small p± , but large values are strongly suppressed by the 
proton form factor, and so the production of states with non-J^ = quantum numbers is 
correspondingly suppressed (see [26, 29] for examples of this in the case of X{c,b) '^{c,h) 
CEP). In particular, we find roughly that 

\n\M = 2)\' {pi? 

|T(J, = 0)|2 ~ (Q2)2 ' ^^-^^ 

which is typically of order ~ 1/50 — 1/100, depending on the central object mass, cms 
energy -y/s and choice of PDF set. As discussed in [17], this 'J^ = selection rule' [32] 
will have important consequences for the case of meson pair CEP. Finally, we note that in 
(2.7) the incoming gluon helicities are averaged over at the amplitude level: this result is in 
complete contrast to a standard inclusive production process where the amplitude squared 
is averaged over all gluon helicities. Eq. (2.7) can be readily generalised to the case of 
non-Jj = gluons which occurs away from the forward proton limit, see in particular 
Section 4.1 (Eq. (41)) of [26], which we make use of throughout to calculate the MM CEP 
amplitude from the corresponding gg — > MM helicity amplitude. 



3. gg vi')vi') amplitudes: Feynman diagram calculation 
3.1 The hard exclusive formalism 

The leading order contributions to the 77 MM process were first calculated in [18] 
(see also [19,33,34]), where M(M) is in this case a flavour non-singlet meson(anti-meson). 
The cross section has been calculated at NLO in [35] . For the case of mesons with flavor- 
singlet Fock states there is also a contribution coming from the LO two-gluon component 
of the meson, and this was calculated in [36] (see also [37,38]). In particular, in [36] they 
considered the process 77 — >■ rjiM (where M ^ rji), where rji is a fiavour-singlct meson 
state that can has both a qq and a gg component. The amplitude can be written as 

M\x'{s,d) = / dxdy(t)i{x)(t)M{y)Tly{x,y;s,e)+ I dxdy (t)G{x)(t)M{y)Tl^,{x,y\s,e) , 
Jo Jo 

(3.1) 

where s is the r/iM invariant mass, x, y are the meson momentum fractions carried by the 
quarks or gluons in the meson. A, A' are the photon helicities and 9 is the scattering angle 
in the 77 cms frame. T^y ^ is the hard scattering amplitude for the parton level process 
77 ~^ <NQQ{99)i see Fig. 2, where each qq and gg pair is coUinear and has the appropriate 
colour, spin, and flavour content projected out to form the parent meson. In the meson 
rest frame, the relative motion of the partons is small: thus for a meson produced with 
large momentum, |/c|, we can neglect the transverse component of the parton momentum, 
q, with respect to /c, and simply write q = xk\n the calculation of T\a'- ^(2^) is the leading 
twist meson distribution amplitude, representing the probability amplitude of finding a 
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(a) (b) 

Figure 2: Representative Feynman diagrams for the 77 — > qqqq{gg) processes. There are 20 
Feynman diagrams of type (a), and the corresponding helicity amplitudes are given in [17] (see also 
[18, 19]). There are 24 diagrams of type (b), and the corresponding helicity amplitudes are given 
by (3.19, 3.20). There are 31 diagrams of type (a), but with the photons replaced by gluons, and 
the corresponding helicity amplitudes are given in [17]. 



valence parton in the meson carrying a longitudinal momentum fraction x of the meson's 
momentum, integrated up to the scale Q over the quark transverse momentum ql (with 
respect to meson momentum k). While (/>Af(x) and <j)i{y) represent the qq distribution 
amplitudes of the mesons M and rji, respectively, (pc corresponds to the gg distribution 
amplitude of the rji. 

We recall (see [17] and references therein for further details) that the meson distribution 
amplitude depends on the (non-perturbative) details of hadronic binding and cannot be 
predicted in perturbation theory. However, the overall normalization of the qq distribution 
amplitude can be set by the meson decay constant /jv/ via [18]^ 

j\xMx) = ^. (3.2) 

It was also shown in [40] that for very large the meson distribution amplitude evolves 
towards the asymptotic form 

4'm{x,Q) -> V3fMx{l-x). (3.3) 

However this logarithmic evolution is very slow and at realistic values the form of cpM 
can in general be quite different. Indeed, the recent BABAR measurement of the pion 
transition form factor Ft^^{Q'^) [41,42], for example, strongly suggests that (f)-,^{x,Q) does 
not have the asymptotic form out to < 40GeV^, although new Belle data [43] are in 
conflict with this. Another possible choice is the 'Chernyak-Zhitnisky' (CZ) form, which 



We note that this normalization, which we use tlirougliout tliis paper, corresponds to the definition 
of /m given in [18] for the case of the pion, that is it with /,r ~ 93 MeV. In the literature the value 
f-n — )■ v^/vr ~ 133 MeV, is often used, in particular in [13,39], which we will refer to later. In Section 5 
of [17], the more conventional value « 133 MeV is confusingly quoted, but in all numerical calculations 
the lower value was correctly taken. 
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we will make use of later on [44] 



^m'{x, = Mo) = sVs/m x{1 - x){2x - if , (3.4) 

where the starting scale is roughly ^uq ~ 1 GeV. For the two-gluon distribution amplitude, 
(poix), the normalization cannot be set as in (3.2), as wc have (pcix) = — (/>g(1 — x), as 
required by the antisymmetry of the pseudoscalar spin projection (3.18) of the gluons under 
this interchange, but an analogous formula can be written down [36, 38, 45] 

f\y4>G{x,Q^)(.2x-l)^fG{Q^). (3.5) 
Jo 

Such an expression serves to define fc, the value of which is to be determined. While 
as ~^ oo, it can be shown that the gg distribution amplitude vanishes due to QCD 
evolution [45,46] 

lim (j)G{x) = , (3.6) 

Q2_).oo 

there is no reason to assume this will be the case at experimentally relevant energies. 

More precisely, the qq flavour-singlet and gg distribution amplitudes can be expanded 
in terms of the Gegenbauer polynomials C„ [40, 45, 46] 

<l>ii,s),Mi^,t^l) = ^77=^(1 - x)[l + E a(^''\„l)C'J\2x - 1)] , 



0g,m(x,/x^) = -4=W^x(1-x) Y1 a^i^^Dct'ii^x-l), (3.7) 

where is the factorization scale, taken as usual to be of the order of the hard scale of 
the process being considered, and Uf = 3 ioi r){') mesons. The f^g (with M = r),i]' in the 
present case) are given by (3.11), with the M dependence expressing the difference due to 

the mixing of the rj, rj states and decay constants. Wc assume that apart from this the 
distribution amplitudes, i.e. the a„, are independent of the meson being considered. We 
choose to include the decay constants ^ fg^ explicitly in our definition of the distribution 
amplitudes, in contrast to, e.g. [13,39], where they are introduced in the hard ampUtude 
Txx' in (3.1)2. 

The evolution of the distribution amplitude is then dictated by the fi^ dependence of 



•^Our choice of normalization of 0i,g differs furtlicr by a factor of l/2\/ Nc, which in [13] is included in 
the hard amplitude, and an additional factor of ^Cf l'infx(\ — x) in (^g due to the different normalization 
of the gg spin projection (3.18). Of course, the final physical result will not depend on the choice of 
convention, which simply corresponds to a choice of which overall factors to include in the distribution 
amplitude ^g,\, and which to include in the hard amplitude r^^' in (3.1). 
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the coefficients a„ via 



7I+V/30 ^,2^ V/3o 

= pi^)ai^^(M^) ( ^ ) + ai-\^^l) ( ^ ) . (3.8) 



That is, the quark and gluon components mix under evolution. Here 
Po = 11 — 2ra//3, the ai^\iJ,Q) at the starting scale /xq are inputs which must be, e.g., 
extracted from data, and the remaining factors 7=*=, p^"^^ are defined in Appendix A. We 
note that despite the different normalization convention we take in (3.7) for the quark and 
gluon distribution amplitudes, the coefficients a„ are defined so that they obey the same 
evolution equation as in for example [13]. While the size of the a^(/i|, = Q'^) will decrease 
with increasing scale, their initial value can in principle be quite large, and the evolution 
of for example the leading coefficient — )• with is slow^. 

Finally, to make contact with the physical rj, rj' states we will be considering in this 
paper, we introduce the flavour-singlet and non-singlet quark basis states 

\qqi) = —^\uu + dd + ss) , 
v3 

\q%) = —i=\uu + dd — 2ss) , (3-9) 
v6 



and the two-gluon state 



\99) , (3.10) 



with corresponding distribution amplitudes given by (3.7). Here, we follow [47] (see also [48, 

49]) in taking a general two-angle mixing scheme for the r] and ij' mesons. That is, the 
mixing of the ij, ij' decay constants is not assumed to follow the usual (one-angle) mixing 
of the states. This is most easily expressed in terms of the rj and r]' decay constants 

fs = fs cos 08 , f1 = -fi sin 01 , 

f^' = fs sin Os , f^' = fi cos 9i , (3.11) 
with the fit of [50] giving 

fs = 1.26/, , 9s = -21.2° , 

/i = 1.17/, , 9i = -9.2° . (3.12) 

We then take the distribution amplitudes (3.7) with the decay constants given as in (3.11), 
for the corresponding Fock components (3.9) and (3.10). That is, the r/ and ij' states are 



^The higher order n = 4, 6,... terms are found to evolve faster towards zero with increasing n. Combined 
with the fact that as n increases, the additional powers of Cl'^ii2x - 1) (or C^/^(2a; - 1) in the qq case) 
give a smaller numerical contribution to the distribution amplitude, this means that any fit can effectively 
truncate the series (3.7) after a limited number of terms. 
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given schematically by (see for example [13, 39] for more details and discussion) 



\v) = /s cos 6*8 



/isineii 4>i,vix,lji')\qqi) + 4>G,r,{x, ixl')\gg) 



\ini') = fssinOs 4>s,n'{x, IJ^f)\<1%) +/icos6'i (t>i,n'{x, liF)\'i'i\) + (t>G,7j'{x, IJ'F)\9a) 



(3.13) 



where to make things explicit the distribution amplitudes (j) are defined as in (3.7), but 
with the decay constants divided out (i.e. (j)s,r]{x, n'jp) = 4)%^j^{x,iJ?p)/ ■■■)■, and these do 
not represent the conventional, normalized expressions for the r/',?] Fock states, but simply 
indicate the decay constants and distribution amplitudes that should be associated with 
each qq and gg state in this two-angle mixing scheme. 



3.2 Preliminary consideration: 77 Mrji 

We begin this section by recalculating the 77 — t- "qiM helicity amplitudes, as in [36]. For 
simplicity, we consider the case of a pure flavour-singlet state ryi, but the following result 
can readilty be applied to the more realistic case of an rj' or rj meson. The appropriate 
qq and gg spin and colour quantum numbers are projected onto the relevant pseudoscalar 
meson states using [39, 40] 



n. = -^=^^, (3.15) 



"c 

where e^^ = — e^^ = 1 and all other components are zero. It can be expressed as 

e^^ = -^e'^'^'^'p^pj , (3.16) 

where e°^^^ = — 1. Explicitly, these correspond to the following combinations of spin states 

1 l' u+{y)v-{l -y) u^{y)v+{l-y) \ _ j5^M f^-.^. 

^ {er{x)e*_^{l -x)- e*_^{x)e:^{l - x)) = -z^ , (3.18) 

where u{z)^v{z) are the usual Dirac spinors for the quarks and e(z) is the polarization 
vector of the gluon, carrying momentum fraction z of the parent mesons, while the it 
signs indicate the particle helicities. The normalizations are conventional (other choices 
are possible, provided the distribution amplitude is suitably re-defined to compensate), 
and we can see that these projections are odd under the parity transformation + ^-t- — , as 
required for a P = — 1 meson state. We make use of (3.14, 3.15) throughout this paper. 
Explicitly calculating the amplitudes corresponding to the 24 contributing Feynman 
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diagrams (see e.g. Fig. 1 (c) of [36])^, we find 



^9(77) ^ j,9{77) ^ , (3.19) 

7.9(77) ^ 7.9(77) I / a^^ - 1 Mir^aelas bcos^ 6 - {2x - l)a 
+- -+ V Nc sy{^-y) a2-62cos2^ ' ^ ' 

where the (77) is to distinguish these from the amplitudes with initial-state gluons that 
we will consider shortly. Cg is the quark charge in the meson M and a, b are given by 

a = {l-x){l-y) + xy , (3.21) 
b={l-x){l-y)-xy . (3.22) 

The '!' in (3.20) indicates that this is not an exact equality, but rather the result of the 
explicit calculation and (3.20) are equivalent after the amplitude T^iff) has been integrated 
over the (antisymmetric) rji distribution amplitude 4>G,r]A^)- That is, they are equivalent 
up to terms which are even under the interchange a; — > 1 — x. We note that (3.19) is 
consistent with the results of [36,37], while (3.20) is consistent up to overall numerical 
factors (in particular our result is a factor of 2 larger than that of [37]). As no definition 
of the gluon spin projector (3.15) is given in [36,37], while the results also differ between 
these two papers by a factor of 2x{l — x), a precise comparison is quite difficult to make, 
and we can reasonably assume that the difference between (3.20) and the results [36,37] is 
due to differing normalization conventions for the spin projectors and gluon wavefunctions. 



3.3 gg ggqq{gg) amplitudes 

Representative diagrams for the three processes we will consider are shown in Fig. 3: in the 
case of the 4:g and 6g amplitudes, we note that all diagrams allowed by colour conservation 
and the antisymmetry of the gluon spin projection contribute to the total amplitude, and 
not just diagrams of this ladder type. After a quite lengthy calculation, we find that the 
gg — > ggqq amplitude for Jz = incoming gluons is given by 



The gg — > gggg amplitude for = incoming gluons is given by 

= T99 = - /^ + ^°ff) (2x - l)(2y - 1) , (3.24) 

We also reproduce for completeness the gg qqqq amplitudes for flavour-singlet mesons 



*We have made use of the FORM symbolic manipulation programme [51] throughout this paper. 



- 10 - 




(a) (b) (c) 

Figure 3: Representative Feynman diagrams for tlic gg — > MM process, wliere the M are flavour- 
singlet mesons. There are 8 Feynman diagrams of type (a), and the corresponding helicity ampli- 
tudes are given by (3.25, 3.26). There are 76 Feynman diagrams of type (b), and the corresponding 
Jz = helicity amplitude is given by (3.23), while a numerical evaluation of the \ Jz\ = 2 case 
is shown in Fig. 4. There are 130 Feynman diagrams of type (c), and the corresponding — Q 
helicity amplitude is given by (3.24), while a numerical evaluation of the | Jj,| =2 case is shown in 
Fig. 4. In the case of the amplitudes for (b) and (c), all diagrams allowed by colour conservation 
are included, and not just diagrams of this ladder type. 



corresponding to the 'ladder diagrams' as in Fig. 3 (a), see [17]. 

^ ^ 64^^«| (1 + cos^g) 

++ — iVc sxy(l-x)(l-y)(l-cos2 0)2 ' ^ > 

rnn- = r'^,- = 647r2a| (l + Scos^g) 

iVc sxy(l-x)(l-2/) 2(1 - cos2 6')2 ■ ^' 

Thus all three = amplitudes are identical, up to overall colour and normalization 
factors (including the factors of '(1 — 2x), (1 — 2y)' which ensure that the amplitudes have 
the correct symmetry under the interchange x(y) -f-)- 1 — x{y)). 

In the case of the gg — )• ggqq and gg — )■ gggg amplitudes for \Jz\ = 2 incoming gluons 
we can find no simple closed form. However, by numerical evaluation we can see that 
they exhibit a similar angular behaviour to the = counterparts. We show this in 
Fig. 4, where we plot the differential cross sections d(TAiA2/d| cos 0| corresponding to the 
amplitudes T^'i_ and T^^_ (or equivalently, and T£^). 

We recall that, due to the selection rule which operates for CEP [32] the contribution 
of the I J^l =2 amplitudes is strongly suppressed. As observed in [17], this is particularly 
important in the case of flavour-non-singlet meson production. These mesons (which 
only have a qq component), cannot be produced in the 'ladder-type' gg — )• qqqq or the 
99 ~^ g9QQ subprocesses discussed above, where the qq pairs forming the mesons come from 
the same quark line. The contributing diagrams are instead of the type shown in Fig. 2 
(left), with the photon pair replaced by gluons, but in this case it was found in [17] that 
the corresponding amplitude vanishes for = incoming gluons, and so the CEP cross 
section for flavour-non-singlet meson production is expected to be strongly suppressed. 

We note that it can readily be shown from the Feynman rules for fermion fields that 
we should associate an additional factor of (—1) with the flavour-non-singlet amplitude of 
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cos^^l Icos^^ 



Figure 4: Differential cross section dcAiAa/dl cos0| for gg — >• 771771 via gg ggqq and gg gggg 
diagrams, shown for illustration, for ^/l = 5 GeV. The qq distribution amplitude is given by (3.4), 
with 71/ = 3, while the gg distribution amplitude is given by (4.4), with /i = f.^. Ai^2 are the 
incoming gluon helicities. 



Fig. 2 (a), for the case of incoming gluons, but not with the 'ladder' diagrams of Fig. 3^. 
That the amplitudes corresponding to these two types of diagrams have a relative minus 
sign is not in fact surprising: if we consider the case that both qq pairs are of the same 
flavour, the diagram of the type shown in Fig. 2 (a) is completely equivalent (replacing 
the photon pair with gluons) to that of Fig. 3 (a) under the interchange of the identical 
fermion quarks (or anti-quarks). 

Finally, it can be shown that the remarkable fact that (3.23), (3.24) and (3.25) are 
identical up to overall colour and normalization factors is not accidental, but can in fact 
be understood in a 'MHV framework. We refer the reader to Appendix B for references 
and a detailed discussion of this. 



4. Results 

From Section 3.3 we can see that to first approximation (ignoring the flavour non-singlet 
component and the \Jz\ =2 amplitudes), the effect of including a gg component to the 77'r7' 
CEP amplitude will simply be to multiply it by an overall normalization factor Ngg, given 
by (omitting the dependence for simplicity) 

^ _ Jdxdy {<Pi,^,{x)cl)i,r,'{y)T'i'{x,y) + 2<Pi^^,{x)My)T'^'ix,y) + M^)<^G^^ 

f dxdyct>i,,,{x)cl),,,,{y)Tii{x,y) 

(4.1) 

with all dependence on the final-state kinematic variables (i.e. cos(0) and s) dropping out, 
due to the identity of the corresponding amplitudes, up to overall colour and normalization 

^In [17], the opposite assignment was taken. However only the quark contribution was considered, in 
which only the relative sign between these two types of diagram was important. The quoted cross sections 
and results are therefore correct. 
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factors. For example for the qq term we may take 

dxdy<^i,,,(x)</)i,,,(y)r^^(x,y) -> -^Uf J dxdy^^0^^0^, (4.2) 

where the factor of (= 3) in the numerator comes from summing over the three quark 
flavours contributing to each qq state, divided by the wavefunction normalizations (3.9). 

To give an initial estimate of the size of this effect, we recall from (3.7) that the gluon 
distribution amplitude can be written in the form 



ct>G,A^,nl) = 0=J^^x{l-x) J2 a^{i,l)C'J_\i2x - 1) . (4.3) 



n=2,4--- 

Neglecting higher order terms in n, we therefore have explicitly 

<f>G,r,'ix,nl) « -^af{fil)x{l - x){2x - 1) . (4.4) 

For the quark distribution, we will consider both the asymptotic (3.3) and CZ (3.4) distri- 
bution amplitudes, i.e. taking a2(/XQ) = and a2{fJ,Q) = 2/3 in (3.7), respectively, with all 
higher n terms being zero. We have 



J x[l-x) 



(4.5) 



d.x^ = 4/l'', (4.6) 

x(l-x) V3 ^ 

r ^ Ja,,,{x){2x-l) ^^j^ , _ 
J x{l-x) 9^6^ 

(4.8) 

Using these we then have, as in (4.2) 

j dxdy<t>rj;^-{x)<i^rr^y)^''^'^^y) 3(/i'')' , 

2 j dxdy<^^7-(x)^G,.'(y)r«^(^,y) ^ ^(/I'? , 

j dxdy<t>'^l,{x)<t^'il,{y)T'i^{x,y) ^ y (/f )2 , 
2 j dxdy<l>^l,{x)ct>GAy)T'''{x,y) ^ ^(/{'? , 
j dxdy^GA^)My)T''{x,y) ^ '^^^{ftf , (4.9) 



where again the n'jp dependence is omitted for simplicity, and it is understood that 0f ^, (x) 
is evaluated at the scale no, so that it is given by (3.4). 
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Figure 5: n = 2 coefficient of tlie expansion (4.3) of the gg distribution amplitude (f>c,M{x, ^"j^), 
as a function of the scale We assume that the flavour-singlet qq distribution amplitude with 
which this mixes under evolution is given by the CZ form (3.4). The dependence exhibited in 
the plot is largely insensitive to the form of the qq distribution amplitude. For < the 
dependence is identical up to the overall factor of (—1). 




To give an initial numerical estimate we may then for illustration take the value 
from [13] 

a^fit(/^o) = 19±5, (4.10) 

which is extracted from the transition form factors F^(/-)^((5^), and where all higher (n = 
4,6...) order terms are neglected. However, this corresponds to the gg distribution am- 
plitude at the scale /Iq = IGeV^. In fact (see also [13]), the leading coefficient a2{fJ-'j?) of 
(/}G{x,fJ-p) evolves rapidly away from its value at fio, decreasing by roughly a factor of ~ 3 
for ^"jp ~ lOGeV^, and decreasing at a much slower rate after that, see Fig. 5. Thus, to 
give an initial numerical estimate we can take CLnit^'p ~ GeV'^) ~ 7, corresponding to an 
experimentally realistic value of Mx ~ 5 GeV. In this case, we find from (4.9) that 

N^/"^ « 9 , (4.11) 
~ 5 , (4.12) 

and therefore we will expect a potentially large increase (~ ^gg) i'^ the Tj(^)rj{') CEP cross 
section for this value (4.10) of a2^(/iQ). While (4.10) may give some rough guidance for the 
expected size of the gg component of the r/('), we note that the fit of [13] contains important 
uncertainties, see the discussion at the end of this section, and so the question of the size 
of (/^o) ^"-"t settled. Therefore, to give a conservative evaluation of the sensitivity of 
the CEP process to the size of this gg component, we will consider in what follows a band 
of cross section predictions, corresponding to the range 

«^ (^0) e (-«?fit/2, +a^fit/2) = (-9.5, 9.5) . (4.13) 

As we will see, even within this quite narrow and conservative range of values, the predicted 
CEP cross section changes considerably. The results we present below are based on an exact 
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numerical calculation, using a modification of the SuperCHIC MC [52], and including the 
precise effect of the evolution of the gg and qq wavefunctions, in contrast to the rough 
estimate described above. We take //q = 1 GeV and = (being of order the 

squared momentum transfer ~ |«|) throughout, although other choices of factorization 
scale are of course possible. 

We show in Figs. 6, 7 the cross sections for ry'ry', -qr] and "qr}' production using the 
CZ form (3.4) for the quark distribution amplitude, for three different choices a2^(/Xo) = 
(—9.5,0,9.5) of the gg distribution amplitude normalization (4.13), and for the mesons 
required to have transverse energy > 2.5 GcV and pscudorapidity |r/| < 1. All cross 
section predictions in this paper are calculated using MSTW08LO PDFs [53]: while there 
is in general a reasonably large variation in the predicted cross sections (which are sensitive 
to the low X and region of the gluon density) on the choice of PDF, we use this set as 
it gives a prediction which is consistent with the existing CDF exclusive 77 data [16], see 
for example [25] for further discussion. We show the cross sections for y/s = 0.9, 1.96 GeV, 
as these are the energies at which the existing CDF data on exclusive final-states have 
been taken, and from a further analysis of which we may hope for an observation of these 
meson pair CEP processes to come [14]. At the LHC, we expect the cross section (for the 
same event selection) to be roughly a factor of ~ 3-5 larger for -y/s = 7-14 TeV, with the 
particle distributions almost unchanged. For o^'CMo) ~ 9.5, we can see that the expected 
cross section increases significantly, while for negative the gq interferes destructively 
with the gg and qq contributions, leading to a suppression in the predicted cross section. 
The predicted 7777 cross section is less suppressed, due to the jJ^j = 2 flavour-non-singlet 
contribution (see the discussion in section 3.3), which for the lower flavour-singlet cross 
sections (corresponding to a2 = 0, —9.5) becomes important. 

We note that the extracted value of a2{fiQ) in general depends on the form of the 
flavour-singlet quark distribution amplitude, i.e. the size of the ai!j(/Xo) in (3.7). With 
this in mind, in Fig. 8 we show the same cross sections at y/s = 1.96 TeV, but for the 
asymptotic form of the quark distribution amplitude (3.3). The ^^'(/^o) ~ ^ cross section 
is predicted to be somewhat smaller than for the CZ distribution amplitude choice, while 
a similar, somewhat larger, relative enhancement is seen for the case that a2'(//o) = 9.5. 
For 02'(/io) = — 9.5 the destructive interference is in fact almost exact for this Mx region, 
and so the cross sections are predicted to be very strongly affected. However, we would 
not necessarily expect to see such a large destructive interference effect in the data, as such 
an exact cancellation may not occur at higher orders. Again, we stress that these only 
correspond to specific cases in a band of possible a2{fiQ), and that the real size of this may 



Neglecting contributions with |Jj| = 2 incoming gluons, the ratio of rj'r]' to rjij {rjVj') 
cross sections are determined by the mixing in (3.11-3.13) to be 



be smaller. 



a{r)'ri') : cr{r]r]') : cr{r]r]) = 1 



2tan2(0i) :tan^(ei) 
1 1 



(4.14) 
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Figure 6: Differential cross section da/dMx for X — rj'jf ,riri,r]T]' production at y/s — 1.96 TeV 
with MSTW08LO PDFs [53], taking the CZ form (3.4) for the quark distribution amphtude, and 
for a band of a^(/i§) values for the gg distribution amplitude. The mesons are required to have 
transverse energy E± > 2.5 GcV and pseudorapidity \ri\ < 1. 



where we have taken the value of 9i from (3.12)^ and the factor of '2' in the 7]7]' case accounts 
for the non-identity of the final-state particles. A measurement of the cross section ratios 
(in particular, a{i]'r]')/a{rjrj')), would therefore serve as a probe of the mixing parameter, 
9i, see (3.12). This ratio is predicted to be unchanged by the inclusion of a non-zero gluon 
component ^ 0, as for = incoming gluons only the flavour-singlet component rji 
of the r] and rj' mesons contributes, while we have seen in Section 3 that the gg — )• ggqq 
and gg — >• gggg amplitudes are identical in form to the purely quark case gg — )• qqqq, and 
will therefore only effect the overall normalization of this flavour-singlet contribution. In 
the ratio of cross sections, this overall factor cancels and we are left with the scaling of 
(4.14) irrespective of the size of the gg component. In Table 1 we show numerical results 
for the ratios (4.14): due to the relative importance of the \Jz\ = 2 flavour non-singlet 
contribution, in the r/r/ case, this scaling is only expected to be approximate, see below 

®We note that the predicted ri'ri and rirj cross sections are lower than in [17]. This is due to the different 
choice of mixing scheme (3.11-3.13) and in particular the lower value of Oi, which leads to a smaller flavour- 
singlet component of the r/. It is found in for example [50,54,55] that this scheme (3.11-3.13) and choice 
of mixing parameters describe the available data well. A measurement of the cross section ratios in (4.14) 
would certainly shed further light on this. 
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Figure 7: Differential cross section da/dMx for X = r]' rj' , -qri , rirj' production at y/s — 0.9 TeV 
with MSTW08LO PDFs [53], taking the CZ form (3.4) for the quark distribution amphtude, and 
for a band of a^(/i§) values for the gg distribution amplitude. The mesons are required to have 
transverse energy E± > 2.5 GeV and pseudorapidity \r]\ < 1. 



(there is also in all cases a small effect due to the differing rj and rj' masses in the MC 
simulation) . 

Thus, to first approximation we can only look at absolute value of the various ??(')??(') 
CEP cross sections to determine the size of the gg component, (Mq). This is potentially 
problematic because of the other uncertainties in the CEP calculation, due primarily to 
the value of the survival factors 5^;^,, S^^^, which are not known precisely, and potential 
higher-order corrections in the hard process, which combined are expected to a give a 
factor of ~ ^2 — 3 uncertainty, as well as a sizeable PDF uncertainty in the low-x, 
regime relevant to such processes (see [25] for a detailed discussion of this and further 
references). Nevertheless, given the sensitivity of the CEP cross section to the gg — t- gggg 
and gg — t- ggqq subprocess, if the gg component of the ?](') is sizeable enough (as the result 
of [13], for example, suggests), such a measurement may still provide useful information. 

However, it is more reliable to look at the ratio of the ??(')??(') cross section to other 
processes, in which case many of the uncertainties due to PDFs and survival factors largely 
cancel out and a potentially much cleaner measurement of the gg component of the r]{') 
becomes possible. With this in mind we show for illustration in Table 2 the ratio of the rjrj 
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Figure 8: Differential cross section da/dMx for X ~ t]' rj' , t]T] , rjrj' production at yG — 1.96 TeV 
with MSTW08LO PDFs [53], taking the asymptotic form (3.3) for the quark distribution amplitude, 
and for a band of a^(/i§) values for the gg distribution amplitude. The mesons are required to have 
transverse energy E± > 2.5 GcV and pscudorapidity \ri\ < 1. 



to tt'^tt'^ and r]'r]' to 77 cross sections'', for the same three choices a2^(/Xg) = (—9.5,0,9.5) 
of the gg distribution amplitude normalization (4.13). 

As mentioned above, the scaling of (4.14) is only approximate, as it ignores the pos- 
sibility of a flavour-non-singlet contribution for the case that the incoming gluons in the 
gg — )• 'r]{')ri{') subprocess are in a | J^l = 2 state. Such a contribution may in particular be 
relevant for rjrj production, when the suppression in the flavour-singlet cross section due to 
the small mixing ~ sin^ 9i can be comparable to the suppression of the flavour-non-singlet 
contribution due to the = selection rule that operates for CEP [32]. This can be seen in 
Table 1, where the ratio of rj'rj' to rjrj (and rjf]') cross sections differs from (4.14) due to this 

^The Tv'^n~ CEP cross section is predicted to be identical to the neutral pion case, up to a factor of 2 due 
to the non-identity of the final-state particles. The observation of such a process has recently been reported 
in [15], although it is not clear that this measurement probes sufficiently high pion k± for the perturbative 
approach discussed here to be applicable. In order to better clarify the situation, a measurement of the 
meson k±- distributions corresponding to the same data would be very useful, while a comparison between 
the k± (and M.^^) distributions at 1.96 TeV and 900 GeV would probe the size of any possible contamination 
due to proton dissociation [3]. We also note that we expect tt+tt" (and K'^K~) data to come soon from 
CMS [56], with a veto applied on any additional particles within their rapidity coverage, which should 
contain a sizeable exclusive component. 
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additional contribution. For larger values of a2'(/^o)i flavour non-singlet contribution 
becomes relatively less important and we approach the expected values; a measurement of 
these ratios may therefore serve as an additional probe of the gg contribution. Recalling 
that the 7r°7r° CEP cross section is also predicted to be strongly suppressed, due to the 
vanishing at LO of the gg — )■ vr^vr" amplitude for = incoming gluons, a measurement 
of the ratios cj(?7(')?7('))/(T(7r°7r°) would also represent as an important probe of the = 
selection rule^. 

We can also see from Table 2 (see also [17]) that the rj'r]' cross section is expected to be 
somewhat larger than for 77 CEP. However, multiplying by the corresponding branching 
ratios, BT{r]' — > 77) 2% and Br(r7 77) ^ 40% [57], squared, we can see in Table 2 that 
the cross sections for r]'r]' (and also for ryry) CEP after branching to the final state are 

predicted to be a small fraction of the direct 77 CEP cross section for the relevant event 
selection, with a similar result holding for the r]r]' final state. Wc therefore do not expect 
these to represent an important background to the existing CDF [16] and any forthcoming 
CMS [58] 77 data^. 
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Table 1: Ratios of vi'M') CEP cross sections at \/i = 1.96 TeV with MSTW08LO PDFs [53], for 

a gg distribution amplitude with difFcrcnt choices of a^(/^Q) and with the qq distribution amplitude 
given by the CZ form (3.4). The meson are required to have transverse energy E± > 2.5 GeV and 
pseudorapidity \t]\ < 1. 
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a(r?r?)/CT(7r'^7rU) 


2.7 


12 


66 




570 


16000 


100000 


CT{ri''n')/a{jj) 


3.5 


100 


660 


a{ri'ri' 47)7(7(77) 


0.0017 


0.049 


0.33 


a{riri 47)7(7(77) 


0.0025 


0.012 


0.066 



Table 2: Ratios of r;(')r/(') to 7r°7r° and 77 CEP cross sections at = 1.96 TcV with MSTW08LO 
PDFs [53], for a gg distribution amplitude with difFcrcnt clioiccs of (i^U-I-q) and with the qq distri- 
bution amplitude given by the CZ form (3.4). The meson/photons arc required to have transverse 
energy E± > 2.5 GeV and pseudorapidity \ri\ < 1. Also show are the ratios <j{ri'r)')/a{'yY) and 
(7(7777)7(7(77) cross sections multiplied by the t]{') 77 branching ratios squared. 



It should be noted that any NNLO corrections or higher twist effects which allow a = contribution 
may cause the precise value of the flavour-non-singlet gg (and tt^tt") cross section to be somewhat larger 
than the leading-order, leading-twist | | = 2 estimate. 

^Moreover there should be a further reduction, in any contamination from g, g' 47 due to the exper- 
imental exclusivity and isolation cuts and event selection. On the other hand it is worth mentioning that 
the multi-photon decay modes of the g and g' are quite sizeable (about 72% and 18% respectively [57]) 
and, in principle, these may also contribute to the background. 
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We recall that while various estimates are available in the literature (see for example [4, 
10,59-62]), no firm consensus exists about the precise size of the gg component^'' of the 
r]{'), which is indicative of the uncertainty present in these 'standard' approaches, see for 
example [10] for a discussion of the theoretical uncertainties ( in e.g. decay form factors) 
present in such extractions. In [13], for example, the gluonic contribution to the r){') 
transition form factor F^(^f^^{Q'^) only enters at NLO. and so has a relatively small effect, 
so that the extracted value relies on a precision fit to the data^^. In the case of CEP, on 
the other hand, the gluonic amplitudes (3.23), (3.24) enter at LO and so, as we have seen 
above, it is potentially strongly sensitive to such a gg contribution. An observation of ryry, 
r/'r/' and/or 7777' CEP could therefore shed important light on this uncertain issue. 

5. Conclusion and Outlook 

In this paper we have performed a detailed analysis of the central exclusive pair produc- 
tion of pseudoscalar rj and rj' mesons, pp ^ p + X -\- p, with X = tjtj, tj'tj' and r/r/'. We 
have concentrated on the case that the mesons are produced at sufficiently high transverse 
momentum k±, that a perturbative approach which combines the Durham model of the 
general CEP process (see e.g. [26] and references therein) with the 'hard exclusive' formal- 
ism (see e.g. [18, 19, 33] and references therein) used to model the specific gg vOvC) 
subprocess, can be applied. 

We have in particular extended the previous work of [17] to include the possibility 
of a non-zero flavour-singlet gg valence component of the 77' (and r/) mesons. We recall 
that knowledge of the quark and gluon components of these states would provide important 
information about various aspects of non-perturbative QCD, see the Introduction for more 
details. With this in mind we have calculated for the first time the amplitudes gg — t- ggqq 
and gg — )■ gggg, where the gg, qq pairs form coUinear pseudoscalar meson states with the 
right spin and colour quantum numbers, see Fig. 3. 

It is worth emphasizing that the parton-level amplitudes relevant to the CEP of 
flavour-singlet (and non-singlet) meson pairs in the perturbative regime exhibit some re- 
markable features. In [17] it was shown that, for the case that the incoming gluons are in a 
Jz = state, the parton-level amplitude for the production of fiavour-non-singlet states, 
99 ~^ QQQQj vanishes completely at LO, while the \Jz\ = 2 amplitude exhibit a 'radiation 
zero' at a particular value of the c.m.s. scattering angle. On the other hand, the production 
of flavour-singlet qq states can take place via a separate set of 'ladder diagrams', see Fig. 3 
(a), and in this case the Jz = amplitude does not vanish, and can be written in a very 
simple form. In this paper, we have shown that if one or both of these qq pairs forming 

^°For example, the situation with regards to the Xc(o,2) decays into r/ and r]' pairs appears somewhat 
puzzling. Experimentally, no enhancement in the production of the r],r]' pairs as compared to pions is 
observed (after taking trivial phase space effects into account). This may indicate that there is some 
destructive interference between the qq and the gg components of the pseudoscalar bosons, or that the gg 
component is small. See [13] for other ideas and [63] for a review and further discussion. 

^'^The potential importance of power corrections to F^(^i^^{Q^) at lower values of Q^, where much of the 
existing data lies, may in fact cast doubt on the reliability of such a fit. We thank Viktor Chernyak for a 
useful discussion on this topic. 



-20- 



the flavour-singlet meson states is replaced by a gluon pair gg, as in e.g. Fig. 3 (b,c), 
then these apparently unrelated amplitudes, gg ggqq and gg gggg, are identical to 
each other and to this 'ladder-type' amplitude, up to overall normalization factors. That 
is, they are predicted to have exactly the same angular dependence in the incoming gg 
rest-frame. In [17] and in Section B we have shown how these remarkable and non-trivial 
results may be explained in the MHV framework by the fact that the same external parton 
orderings contribute in all three cases. 

In CEP, this contribution from a gg valence component of the r/', 7] mesons enters at 
the same (leading) order to the purely qq contribution. We have seen in Section 4 that 
as a result the predicted ??(' )??(') CEP cross sections display quite a strong sensitivity to 
any potential gg component of the rj' , r] mesons. An observation of r/ry, rj'rj' and/or rjij' 
CEP could therefore provide a potentially powerful handle on the structure of the r]', r] 
states. More generally, the observation of such meson pair {rjOrjC), vrvr, KK...) CEP 
processes would serve as a test of the unique properties of the gg — )• MM parton-level 
helicity amplitudes, in particular the vanishing of the = amplitudes for flavour non- 
singlet mesons and the identical form of the gg and qq mediated flavour-singlet amplitudes. 
We emphasize in particular that the ij'r]' CEP cross section is predicted to be very large, 
see Section 4. These results, in conjunction with the = CEP selection rule, lead to 
highly non-trivial predictions for the meson pair CEP cross sections, which it would be 
very interesting to compare to data. 

Finally, we note that we may expect new results on ri{')r]{') and meson pair CEP 
more generally to come from further analysis of the existing CDF data (in particular the 
existing 47 candidates with Et > 2.5 GeV and forward rapidity gaps) as well as from 
the CMS/Totem (ATLAS) special low-pileup runs with sufficient luminosity [14, 15]. At 
CDF, the observation of 77 CEP has already been reported [16], and it was determined 
experimentally that the contamination caused by tt^tt^ 47 CEP, with the photons in the 
7r° decay merging or one photon being undetected, is very small (< 15 events, corresponding 
to a ratio N{TT^Tr'^)/N{^^) < 0.35, at 95% C.L.), in agreement with the prediction of [17], 
for which a{'K^'7T^)/a{'y'y) ~ 1%. This prediction is a non-trivial result of the perturbative 
CEP framework and the hard exclusive formalism, and an observation of vr'^Tr'^ CEP, which 
may hopefully come with the increased statistics that a further analysis of the existing 
data can bring [14], would certainly represent an interesting further test of the theoretical 
formalism. In Section 4 we observed that the rj'r]' and rjr] to 47 cross sections are expected 
to be of about the same size or maybe even be larger than the tt^tt*^ cross section, but 
that they should not be an important background to 77 CEP. Nonetheless, this raises the 
possibility of a future observation of these processes via the 77 decay chain. We also note 
that the branching ratios for the p^j and -qn^TT^ decays of the 7]' are sizeable, and may 
be viable channels for the observation of rj'ri' (and ryry') CEP, a possibility which would 
hopefully be confirmed after further analysis and simulations. 

To conclude, the CEP of r]' and 77 meson pairs, in the perturbative regime, represents 
a novel (and complementary) probe of the size of a flavour-singlet gg component of these 
mesons. It is our hope that future rj' , 77 pair CEP data and analysis will be forthcoming 
from the Tevatron and the LHC and that through this we can shed some light on this 
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interesting and currently uncertain question. 
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A. Meson distribution amplitudes: anomalous dimensions 

The anomalous dimensions which control the evolution (3.8) are given by [45,46] 
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In terms of these we can then write down the eigenvalues 7^^^ which diagonalise the anoma- 
lous dimensions 
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The parameters in (3.8) are given by 
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B. MHV calculation 



It is well known (see for example [64]) that the tree level n-gluon scattering amplitudes, in 
which the maximal number (n— 2) of gluons have the same helicity, the so-called 'maximally 
helicity violating' (MHV), or 'Parke-Taylor', amplitudes, are given by remarkably simple 
formulae [65,66]. These results were extended using supersymmetric Ward identities to 
include amplitudes with one and two quark-antiquark pairs in [64], where 'MHV refers 
to the case where (n — 2) partons have the same helicity. In these cases, simple analytic 
expressions can again be written down for the MHV amplitudes, while for greater than 
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2 fermion-anti-fcrmion pairs (recalling that the helicities of a connected fermion— anti— 
fermion pair must be opposite) no MHV amplitudes exist. 

With this in mind, we will show that the fact that (3.23), (3.24) and (3.25) are iden- 
tical up to overall colour and normalization factors is not accidental, but follows from the 
observation that the Jz = helicity amplitudes considered in Section 3.3 are MHV, with 
n — 2 = 4 partons (the two incoming gluons, and two outgoing partons) having the same 
helicity. We will in particular show that the same set of MHV partial amplitudes, i.e. with 
the same orderings of the parton momenta, contributes in all cases. This section explores 
purely theoretical aspects of the previously calculated amplitudes, and so a reader who is 
only interested in the CEP cross section predictions may skip forward to Section 4. 

We recall that in general it is well known that a full n-parton amplitude Ain can be 
written in the form of a 'dual expansion', as a sum of products of colour factors r„ and 
purely kinematic partial amplitudes An 

Mn{{Pi,hi,Ci}) = ig"-^J2'^n{(r{ci})An{a{l^\--- ,n^") , (B.l) 

cr 

where q are colour labels, i^* corresponds to the ith particle {i = 1- ■ ■ n), with momentum 
Pi and helicity Aj, and the sum is over appropriate simultaneous non-cyclic permutations 
a of colour labels and kinematics variables. The purely kinematic part of the amplitude 
An encodes all the non-trivial information about the full amplitude, Mn, while the factors 
r„ are given by known colour traces, see for instance [67] for more details. Adjusting the 
notation from (B.l) slightly for clarity, the n-gluon and qq {{n — 2)-gluon) MHV partial 
amplitudes are given by 

Aig+,g+,...,g-,...,qJ,q+^^,...,g+) = ^n^_^^(^jjj.\^l^ ' (^-3) 

where {ki kj) = {k'~\k'^) = u-{ki)u+{kj) = v+{ki)v-{kj) is the standard spinor contraction, 
and all momenta are defined as incoming. For the case that the quark (anti-quark) has 
positive (negative) helicity, it is sufficient to simply interchange j and j' + l in the numerator 
of the right hand side. The form of (B.3) expresses the important requirement that gluons 
are always emitted from the 'same side' of the connected quark-antiquark line (see for 
example Fig. 1 of [68] and the discussion in the text); that is, the qq pair must appear 
consecutively and in the same order. 

We can then apply these general formulae to the specific n = 6 amplitudes in Sec- 
tion 3.3, where the outgoing qq and gg pairs must form coUinear meson states of the 
correct colour and spin. Firstly, we can see that in the case of both of the n parton am- 
plitudes (B.2) and (B.3), the numerators arc independent of the particular (non-cyclic) 
permutation of the partons being considered, and will therefore factorize. The form this 
takes is given by the corresponding spin projections, as in (3.17) and (3.18). In particular, 
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if we define the following momenta 



^3 = xp3 Z4 = (1 - x)p3 I5 = ypi 1% = {1- y)pA , (B.4) 

where 393^4 are the 4-momenta of the outgoing mesons, then the numerator in the 6-gluon 
case corresponds to (recalling that the helicity here is defined with respect to the incoming 
gluon momenta) 

^(4 6)^ + (3 5)^ -(3 6)^ -(4 5)^ 

= |2(2y-l)(2x-l) . (B.5) 

In the qq case, the final expression is the same, but with the factor of '(2?/ — 1)' removed. 

We can see from (B.2) and (B.3) that the amplitudes for a given ordering of partons 
are identical between the two cases, up to these overall numerator factors. To justify the 
statement that it is indeed the same set of orderings which contribute in the case of both 
amplitudes (3.23) and (3.24), we must consider the colour factors, T„, in (B.l). In the qq 
case this is given by 

r„((n - 2)g + qq) = (A^ • • • X'^-^.j, , (B.6) 
where ii{ji) is the colour index the quark (anti-quark). In the ra-gluon case, T„ is given 

by 

r„(n5) = Tr(Ai---A«) . (B.7) 

For a given T„, the corresponding diagram has the same cyclic ordering of the quark and 
gluons as their colour labels in T'^. The A matrices are normalized (unconventionally) so 
that Tr(A"A*) = 5"'', as required by the form of the expansion given in (B.l). 

Now, an explicit calculation in Section 3.3 has shown that in both cases it is only the 
leading terms in Nq which give a non-zero contribution. In the gg — > qqgg case the terms 
of order 1/Nc which come from certain 'abelian diagrams' (i.e. those which are identical 
to the 77 — >■ qqgg diagrams, with the photons replaced by gluons) do not contribute due to 
the vanishing of the = two-photon amplitudes (3.19), while in the gg — t- 4:g case, all 
colour factors coming from the colour-singlet final and initial states are of order N^. We 
therefore know that all the partial amplitudes with colour factors T„ which are sub-leading 
in Nc must sum to give zero contribution, and can be neglected. By considering the various 
possible values of for different non-cyclic orderings of the quarks and gluons, it is easy 
to show that only those diagrams of the type shown in Fig. 9 (d-f) and Fig. 9 (g— i) give 
such a leading Nc colour factor, with in particular 

T6{4g + qq) = Nc + --- , (B.8) 
n{6g) =Nl + -- - , (B.9) 

in all cases. That is, the same type of diagrams (corresponding to a given ordering of the 
external parton momenta) contribute in both cases. Recalling that the denominator term 
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Figure 9: Representative contributing MHV diagrams for flavour singlet meson pair production, 
via (a-c) gg qqqq (d-f) gg qqgg (g-i) gg gggg 



in (B.2) and (B.3) is identical for a given particle ordering, while the numerator and colour 
factors factorize, we can see that the resulting amplitudes will be equivalent up to these 
overall factors, as expected. 

Explicitly, calculating the partial amplitudes corresponding to the 6g diagrams shown 
in Fig. 9 (g-i), including those coming from interchanging pi ■<->■ p2, x — >■ 1—x and y — >■ 1—y, 
we find 

At' = -2^^^ , (B.IO) 



su 
s 



2*2 



J^^ = , (B.12) 

where these correspond to the kinematic amplitudes with the numerator factors, as in 
(B.5), as well as an additional factor of xy{l — x)(l — y) in the denominator, omitted for 
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simplicity, and the subscript ('g,h,i') indicates the corresponding diagram in Fig. 9. The 
amphtudes in the Ag case are all a factor of 2 smaller due to the requirement that only 
those diagrams with the gluons emitted on one side of the quark line are included, while 
there is no such requirement for the purely gluonic case. 

In fact, the individual MHV amplitudes corresponding to these parton orderings are 
divergent, due to factors of {hh) ~ = and {hie) ~ = present in the denominators, 
but the sum over all particle interchanges (in particular, of the two gluons forming the 
meson states, i.e. hih) O hih), where there is only one pair to be interchanged in the 
ggqq case) is not. To arrive at these (finite) results requires a careful grouping of the 
contributing amplitudes and application of the Schouten identity 

(12) (3 4) = (14) (3 2) + (13) (2 4) . (B.13) 

Combining (B.12) with the relevant numerator, colour and normalization factors, we can 
readily reproduce the expressions of (3.24) and (3.23), up to an overall phase, which we 
take to be due the differing convention for the overall phase of the amplitude which the 
normalization of (B.l) implies. 

Finally we turn to the gg — )■ qqqq amplitude. We first consider the colour factor, T„, 
in (B.l) for the qqqq case. This is given by 

r„((n - 4)g + 2qq) = ^^(A"^ • • • A"Onai(A'^ • • • X'^i^a, , (B.14) 

where ii,i2 are the colour indices of the quarks, ai,a2 are the colour indices of the anti- 
quarks and the labels a.j, 5j refer to the gluons. The pair (0102) is a permutation of the 
pair (J1J2)) where quark is connected by a fermion line to antiquark j^. p is then the 
number of times ak = jk, with p = 1 if (aiQ;2) = (iij2)- The sum is over all the partitions 
of the gluons {I + I' = n — 4, Z = 0, • • • , n — 4) and over the permutations of the gluon 
indices, with the product of zero A matrices given by a Kronecker delta. 

In the case of (3.25), we can again see that it is only the leading terms in Nc (0(1)), 
which contribute to the hard amplitude. Again, it can readily be shown that the relevant 
partial amplitudes for this have the same parton ordering as those in the 6g and 4g cases, 
as shown in Fig. 9 (a-c). Diagrams a and b correspond to the p = and c to p = 1 in 
(B.14), with the additional factor of 1/Nc in the case of diagram c canceled by the factor of 
6ii in the numerator; in both cases, the colour factor is then 0(1). In the case of diagram c 
both gluons must be emitted by one of the quark lines, otherwise the amplitude will vanish 
due a factor of Tr(A) in the numerator. 

The kinematic amplitude for two non-identical quark anti-quark pairs, qq and QQ, is 
given by (see e.g. [69]) 
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where 



= {mnf 
F(-,+) = (lm)2 



F(+,-) = (nm + l)2, 
= {lm + lf . 



(B.16) 



We note that (B.15) and (B.16) are for the case that p = in (B.14), i.e. for diagrams a 
and b of Fig. 9. For the {p = 1) case of diagram c, we must simply interchange n m 
in the above expressions. We can again readily show that summing over the contributing 
quark helicity states according to the spin projection (3.17), as was done in (B.5) for the 
65 case, will simply give an overall factor of for the numerator of (B.15), irrespective of 
the particle ordering. 

Thus, the denominator term in (B.15) is identical to that in (B.2) and (B.3) for a 
given particle ordering, while the numerator and colour factors factorize. As it is again the 
same type of diagrams (corresponding to a given ordering of the external parton momenta) 
which contributes, we find that the resulting amplitudes will be equivalent up to these 
overall factors, as expected^^. 
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